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Work distribution function for a Brownian particle driven by a nonconservative force 
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We derive the distribution function of work performed by a harmonic force acting on a uniformly 
dragged Brownian particle subjected to a rotational torque. Following the Onsager and Machlup’s 
functional integral approach, we obtain the transition probability of finding the Brownian particle at 
a particular position at time t given that it started the journey from a specific location at an earlier 
time. The difference between the forward and the time-reversed form of the generalized Onsager- 
Machlup’s Lagrangian is identified as the rate of medium entropy production which further helps us 
develop the stochastic thermodynamics formalism for our model. The probability distribution for 
the work done by the harmonic trap is evaluated for an equilibrium initial condition. Although this 
distribution has a Gaussian form, it is found that the distribution does not satisfy the conventional 
work fluctuation theorem. 


1. INTRODUCTION 


Fluctuations are ubiquitous in most of the physical 
processes we observe in nature. In general, though the 
effect of fluctuations could be ignored for the coarse 
grained description of the average thermodynamic quan¬ 
tities, these are no longer negligible for small systems, 
like, biomolecules, nano-machines etc. For instance, the 
movement of a molecular motor of nanometer size is af¬ 
fected signihcantly by the mean thermal energy of the 
surrounding medium which causes the motor’s position 
to fluctuate. In linear irreversible processes, these fluctu¬ 
ations obey a general principle known as the fluctuation- 
dissipation theorem [1, 2]. However, in the case of sys¬ 
tems far away from equilibrium, the assumptions of the 
linear-response theory are not valid since the fluctuations 
of the thermodynamic variables do not follow the linear 
law of relaxation. 

Since the finding of the fluctuation relations for the 
probabilities of observing the entropy producing and con¬ 
suming trajectories in simulations of shearing fluid nearly 
two decades ago [3], there has been a vast number of 
studies [4-19] extending the ideas and formalisms of lin¬ 
ear irreversible thermodynamics to nonequilibrium sys¬ 
tems. The emphasis of the fluctuation relations has 
been on quantifying the asymmetry between the prob¬ 
ability distributions of a time-integrated quantity St to 
have values ±IF in the forward and time-reversed pro¬ 
cesses, and it is expressed in the mathematical form as 
= exp(VF/fcsT). St could represent the 
work done by an external force, dissipated heat, entropy 
production etc, and sj stands for the same thermody¬ 
namic quantity in the time-reversed process. These re¬ 
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lations, valid for systems far-from-equilibrium, could be 
reduced to the fluctuation-dissipation theorem for near¬ 
equilibrium processes [7]. Of the many ramifications of 
the fluctuation relations, the most notable results are the 
Jarzynski equality [8] and Crooks work relation [9] which 
relate the asymmetry to the free energy difference be¬ 
tween the initial and hnal equilibrium states which are 
connected through a reversible or irreversible transfor¬ 
mation. Although, to prove the fluctuation relations, one 
does not always require to obtain the detailed distribu¬ 
tion functions, to verify them experimentally, however, it 
is preferable to get the histograms of the relevant thermo¬ 
dynamic quantities. In many situations [20], a detailed 
knowledge of the distribution of the nonequilibrium work 
might be essential to compute the equilibrium free en¬ 
ergy difference utilizing the Jarzynski equality. Despite 
the significant progress in the theoretical studies of these 
distribution functions for work, heat, entropy produc¬ 
tion etc for a variety of nonequilibrium systems [17-29] 
(and the references therein), the derivation of the exact 
analytical expressions for these distributions are not so 
many. This is possibly due to the increasing difficulty in 
obtaining the analytical solutions in closed forms as the 
dynamics of the system becomes more complex. 

In this paper, we derive the positional probability dis¬ 
tribution and the work distribution function for a Brow¬ 
nian particle using the Onsager-Machlup’s functional in¬ 
tegral approach. In their theory of linear relaxation pro¬ 
cesses, Onsager and Machlup developed a path integral 
approach [30, 31] for the probability distribution func¬ 
tions of trajectories in the configuration space with spe¬ 
cific initial and hnal states. Using a variational technique 
they subsequently showed that this approach leads to On- 
sager’s minimal energy dissipation. Recently, there have 
been many efforts to generalize the Onsager-Machlup 
theory for systems in nonequilibrium steady states [32- 
39]. Extending this approach, Taniguchi and Cohen 
[36, 37] obtained the distribution functions of work, dis¬ 
sipated heat for a uniformly dragged Brownian parti¬ 
cle. Further, after expressing the Onsager-Machlup La- 
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grangian in terms of the dissipation functions and the 
entropy production rate, they showed how these fluctu¬ 
ating quantities satisfy the energy conservation law and 
the second law of thermodynamics. Similar variational 
scheme has also been considered to obtain the probability 
distributions of thermodynamic quantities for a variety 
of systems [40]. 

In the present work, we consider two dimensional mo¬ 
tion of a Brownian particle which is subjected to a non¬ 
conservative force and a confining harmonic potential 
with its minimum moving with a uniform velocity v. The 
nonconservative force 

/ = -ynx + xny, (1) 

with n denoting the strength of the force, produces an 
anticlockwise drift about the origin. The nonconserva¬ 
tive force can arise due to a vortex centered at the ori¬ 
gin of the system in the surrounding liquid [41]. Both 
the time dependent confining potential and the noncon¬ 
servative force are responsible for breaking the detailed 
balance and driving the system out of equilibrium. For 
this model, we derive the joint probability distribution 
which gives the probability of finding a Brownian parti¬ 
cle at a given position at a specific time. After defining 
the forward and the time-reversed trajectories, we ob¬ 
tain the expression for the dissipated heat and the en¬ 
ergy conservation relation. In the process, we identify 
various works, namely, the work done by the harmonic 
trap which can be identified as the thermodynamic work 
[8] and the mechanical work done by the nonconserva¬ 
tive force [28, 29]. Finally, using a variational approach, 
we obtain the distribution function for the work done by 
the trap. The work distribution has a Gaussian form 
with the mean and the variance having a nontrivial oscil¬ 
latory dependence originating from the nonconservative 
force. Although, this distribution is Gaussian, it does not 
satisfy the conventional fluctuation relation. Physically, 
the presence of the nonconservative force makes the role 
of the thermal fluctuation negligible. Similar deviation 
can also be found in [42, 43] for the motion of a torsion 
pendulum driven by an external time dependent torque. 

The rest of the paper is organized in the following way. 
In section 2, we introduce two coupled Langevin equa¬ 
tions describing the motion of the Brownian particle in 
the X — y plane. In section 3, we hnd out the joint prob¬ 
ability distribution for the position of the Brownian par¬ 
ticle. Starting with the Fokker-Planck equation for the 
joint probability distribution, we express first the transi¬ 
tion probability in a functional integral form and subse¬ 
quently evaluate the path integrals employing Onsager- 
Machlup’s variational approach. The joint probability 
distribution is obtained from the expression of the transi¬ 
tion probability. From the generalized Onsager-Machlup 
Lagrangian, we define in section 4 the thermodynamic 
work, mechanical work, dissipated heat and the entropy 
production of the medium. We further show here the va¬ 
lidity of the energy conservation relation and the second 
law of thermodynamics for a trajectory. In section 5, we 


obtain the explicit form of the distribution function of 
the thermodynamic work. We summarize our findings 
and discuss possible future research on this model in sec¬ 
tion 6. Some of the details of the derivations and algebra 
are presented in the appendices. 


2. THE MODEL 

The Brownian particle considered here moves 
in a plane under the influence of the potential 
U[x{t),y{t),t] = I {(a;(t) — vt)"^ + y{t)‘^} and a noncon¬ 
servative force as given in equation (1) [29, 41]. Here v 
is the constant velocity with which the harmonic trap 
moves along the x direction. 

The overdamped Langevin equations describing the 
motion of the particle in the x-y plane are expressed as 

X = —^-{x - vt) - ay+ (2) 

Tr a 

y = -—y + crx + -^y{t), (3) 

T'-p Ct 

where, a is the viscous drag coefficient, a = n/a and 
Tr = a/k is the relaxation time of the trap. The overdots 
imply derivatives with respect to time, t. ^x{t) and ^y(t) 
represent the Gaussian-distributed thermal noise having 
zero mean and correlations as 

2ry 

= ( 4 ) 

Here (3 = {kBT)~^ and f, j correspond to the x, y coor¬ 
dinates. 


3. PROBABILITY DISTRIBUTION 
FUNCTIONS RELATED TO THE LOCATION OF 
THE PARTICLE 


The joint probability distribution is defined as the 
probability of finding the Brownian particle at position 
{x, y) at time t. The time evolution of the joint probabil¬ 
ity distribution is governed by the Fokker-Planck equa¬ 
tion which can be written as 


dp{x,y,t) 

dt 


^FP p{x,y,t), 


( 5 ) 


where, the Fokker-Planck operator Lpp is given by, 


d d 

Lpp = -^{{x-vt)/Tr + ay} - —{-y/Tr + ax} 


+D 


dx"^ 


D 




( 6 ) 


with the diffusion constant, D = l/(a/3). Let us con¬ 
sider the simplest situation where v = a = Q. Writing 
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the Fokker-Planck equation in terms of the probability 
current, Jp, we have 


dt 


-V-Jp, 


( 7 ) 


where 


Jp — X 


xp Ddp 


yp Ddp 

Tr dx 

"T y 

Tr dy 


( 8 ) 


Substituting p = exp[—/3C/] (with v = 0), we obtain a 
vanishing stationary current corresponding to the equi¬ 
librium situation of a fixed harmonic potential. This is 
the situation where the detailed balance is satisfied. The 
nonconservative force which cannot be expressed in the 
form of a gradient of a potential violates the detailed 
balance. 

p{x, y, t) can be expressed in terms of the transition 
probability P{x,y,t \ XQ,yo,tQ) which gives the proba¬ 
bility of finding a Brownian particle at position {x, y) at 
time t, given that it was at (a:o,yo) at an initial time to- 
For simplicity, we choose to = 0 and denote the transi¬ 
tion probability as P{x,y,t \ Xo,yo). The Fokker-Planck 
equation (5) allows us to obtain the following path in¬ 
tegral representation for the transition probability [44] 
(some of the details of the derivation are presented in 
appendix- A) 


P{xf,yf,tf I xo,yo) = j Pixit)] j V[y{t)] 


X exp 


pf 


P{x,y,x,y,t) dt 


(9) 


C{x, y, X, y, t) is the generalized Onsager-Machlup La- 
grangian having the form 

^x,y,x,y,t) = -^[{x + {x - vt)/Tr + ayf 

+ {y + y/rr-crxf]. ( 10 ) 

We evaluate the functional integrals (9) using a varia¬ 
tional approach [30, 31, 36, 37, 45], in which one evaluates 
the most probable path that gives the most significant 
contribution to the functional integral. This is expected 
to be a good approximation for low noise strength or 
H —>■ 0 [45]. The approximate form of the transition 
probability can be written as 


I xo,yo) oc exp(- 




x,y,x,y,t) dt 


( 11 ) 


where the object in the square bracket is evaluated over 
the most probable path leading to the largest value of 
the exponential. The most-probable path is determined 
by using the extremization condition 

s[ jC{x,y,x,y,t)dt = 0. ( 12 ) 

^0 


Doing a Taylor-expansion around the optimal values, 
{xc,yc), and retaining fluctuations about the optimal 
path, u = (x — Xc) and z = {y — yc) up to second or¬ 
der the transition probability in (9) can be written as, 

P{xf,yf,tf I xo,yo) = 

xexp^-J C{xc,yc,Xc,yc,t)dt^ . 

(13) 

Since we consider trajectories with fixed end-points and 
allow no variation at the end points, u and z sat¬ 
isfy Dirichlet boundary conditions u(0) = 0, u{tf) = 
0, z(0) = 0, z{tf) =0. As a consequence of this, the 
fluctuation factor 

= / PW)] / P[z{t)] 

Ju{0)=0 Jz{0)=0 

X ([7* Pu/Tr azf' + [z + z/Tt — crw]^) |, 

(14) 


is a function of the final time t f only. Instead of evaluat¬ 
ing the path integrals in Eq. (14) directly, we determine 
4>{tf) using the normalization condition of the transition 
probability, i.e., the total probability of finding the Brow¬ 
nian particle at any final value of (a:/, t//) after time t/ is 
unity. It can be proved, a posteriori, from the normaliza¬ 
tion property of P{xf,yf,tf \ XQ,yQ) |„=oj that the final 
expression of (j){tf) has the form 

The extremization condition (12) leads to two coupled 
Euler-Lagrange equations for the optimal paths Xc{t) and 
yc{t) 

Xc - (ct^ + \ItI)xc P 2ayc Px[\ - — ) = 0 (16) 

V A Xr J 

jjc - [cP' -I- \/T‘^)yc - 2aXc P vat/xr = 0. (17) 

The solutions of the above equations are 

Xc{t) = v[bt — Tr P a^T^]/P 

P cosh(t/T,r)[a 2 cos(crt) — fli sin(crt)] 

-I-sinh(t/rr) [04 cos(crt) — a 3 sin((Tt)], (18) 


yc{t) = CFTrV[bt — 2 Tr]/P 

P cosh(t/Ti.)[ai cos(crt) -I- 02 sin(crt)] 
-|-sinh(t/Tr.)[a 3 cos(crt) + 04 sin((Tt)], (19) 

where 6 = (1 -I- cr'^r^)- The constants Oi, 02 , 03 
and 04 are calculated using the boundary conditions 
a;c(0) = xo, Xc{tf) = xj, yc{d) = yo, ydtf) = yj. These 
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constants are To obtain the joint probability distribution 

p{xf,yf,tf), we integrate the above expression over the 
= yo +‘2,crvT^ /b , (20) initial points, assuming that the system was initially in 

02 = Xq — v{a^T^ — Tr)/b'^, (21) thermal equilibrium with the joint probability distribu- 

^3 = -^^^^^^^|cos(CTt/) (-2ar2u + 6t/wr,-5V) p{xo,yo) = § exp ^ j/g]). So, 

+ (2crT,^u + fc^yo) cosh(t//rr) 

+ sm{atf) (^TrV — a^T^v — btfv + j^, (22) 


04 = 


CSCh(t//Tr-) 


^2 ■^cos{atf){TrV-a^T^v-btfV + b^Xf) p{xf,yf,tf) 

+ COSh(t//rr) { — TrV + O^T^V — b^X q) 

+ sin((Tt/) {2(jt^v — uTrbtfv + 6^y/) j-. 


= dxQ dyopixo,yo)P{xf,yf,tf \ xo,yo) 


(23) 


27r 


exp 


/?/ce 


-2tf/Tr. 


262 


X (^v'^T^ + + 2e^fP^vTrqj , (25) 


Substituting Eqs.(14), (18) and (19) into Eq.(13), and 
performing the integration over time, we have the com¬ 
plete expressions for P{xf,yf,tf \ xo,yo) as 


X exp 


I xo,yo) = 

(- 1 - 1 - 


[1 -I- coth(f//Tr)] 


SDXr 


dnDTr 

[(oi -|- 03)^ -|- (02 -l- 04)^] [■•( 24 ) 


where 


0 . 4 ^ 



FIG. 1: Plot of the joint distribution function 
p{xf,yf,tf) for the parameter values D = 0.05, /3 = 10, 
Tr = 10, k = 0.2, and z; = 2 at final time tf = 20. 


p = v'^t'jb + 2TrVXf — 2a‘^T^vxf + {x'j + y‘j){l + a^T^) 
+t'^[v^ + davyf + 2a'^{x'j + yj)] - 2tfVXf 
— 2tfV^Tr — 2tfVTra[yf aTr{xf -I- CTTf.?//)], (26) 

q = cos{atf)[vtf — Xf + Tr{—v + a^Tj-Xf — 2ayf)] 

+ sm{atf)[—yf + <7Tr{—vtf + 2xf + aTry/)]. (27) 


During simplifications, we have used the relation, 
PDkxr = 1. For a particular set of parameter values, the 
joint probability density p{xf, yf,tf) is plotted in Fig.(l). 

The marginal probability is defined as the probabil¬ 
ity of finding the particle at a given Xf or yf for arbi¬ 
trary values of yf or x/, respectively. After integrating 
p{xf,yf,tf) with respect to the appropriate end point 
variable, the marginal probability distributions are found 
as 
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P{xf,tf) = I dyf p{xf,yf,tf) 
27r 


'V'" / ' 

-I exp| — 

X ti — Tr + cr^T^] + + fTr [(—1 + cr^r^) cos{crtf) + 2aTr sin(cri/)]^ (28) 


and 


PiVf^tf) = / dx/ p{xf,yf,tf) 




1/2 


exp< - 


kp e 




271 y ( 2&4 

X wcrr,. [t/6 — 2rr] + y/ 6^ + iiTr [—2(Trr cos(cri/) + (—1 + cr^ry) sin((Ti/)]^ (29) 

I- 


Eqs. (25), (28) and (29) are the main results of this 4.1. Entropy production in the medium 

section. 


4. GENERALIZED ONSAGER-MACHLUP 
LAGRANGIAN AND STOCHASTIC 
THERMODYNAMIC QUANTITIES 

In this section, we establish an energy-balance rela¬ 
tion and a second law type formulation valid for a sin¬ 
gle trajectory. We consider that, for every forward tra¬ 
jectory [x{t),y{t)] specifying the system’s evolution over 
time interval [—t, r], there exists a time-reversed trajec¬ 
tory [a;f (t), yf (t)] over the same interval. The probability 
measure of the forward path could be expressed in terms 
of the transition probability up to a normalization factor 
as. 


The ratio of the transition probability for the forward 
trajectory to that of the reversed trajectory is given by 


PF(x(T),y(T) I x(-T),y(-r)) 
PR{x^T),y^T) I a:t(-T),yt(-T)) 


= exp 


-J {C-C'‘)dt 


(32) 


where 

L — = P [kx{x — vt) -I- kyy + nxy — nyx] . (33) 

The entropy production in the medium ASm satisfies the 
detailed fluctuation theorem [28] 


PF{x{T),y(T) I x{-T),y{-T)) 
^exp(^-J C{x,y,x,y,t)dt^ , 


(30) 


with the generalized Onsager-Machlup Lagrangian 
C{x,y,x,y,t) as given in Eq.(lO). Similarly, the transi¬ 
tion probability corresponding to the time reversed tra¬ 
jectory is expressed as [36, 37] 


Pr{x'^ iT),y'^ (t) 1 x\-T),y\-T)) 


exp 


C\x\y\x\y\t) dt 


(31) 


AS„ 


, , / PFix{T),y{T) I x(-T),y(-T)) 

^ \PR{x^T),y^T) \ x'f{-T),y'^{-T)) 


From relation (32) and (33), we have 





T 


-AU + 




dt -b y {—nxy + nyx)dt , 
(35) 


where AU = U{x{T),y{T)) — 17(a:(—r), y(—r)) is the 
change in the potential energy and the explicit form of 
the second term is 


where is obtained from C under the time-reversal 
changes v —>■ —v, x^t) —>• x{—t), x'^{t) —>• 
-x{-t), y^t) y{-t) and y\t) -y{-t). 



dt 


dt = —kv 



— vt). 


(36) 
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4.2. Work and dissipated heat 

Since the heat transferred to the medium is Q = 
TASm, we use Eq. (35) to find 

r du r 

—AU + J dt + J {—nxy + nyx) dt 


Q = 

= -ACZ + Wtr + Wm. 


(37) 


In the last equality, the dissipated heat is expressed in 
terms of the change in the potential energy and the work 
that consists of two parts, namely, the thermodynamic 
work, Wtr = ^ dt, which is the work done on the 

system by the harmonic trap [18, 19, 21, 46] and the me¬ 
chanical work, Wni = ^^^{—nxy -I- nyx) dt. Similar par¬ 
tition of the total work into physically meaningful works 
can be found in recent literatures [42, 47, 48]. For further 
discussions on various definitions of work and their be¬ 
havior under gauge transformation, readers are referred 
to [48]. Together with these definitions, Eq.(37) becomes 
the statement of the energy conservation principle valid 
on each trajectory [13, 14, 46]. 

Next, to show the second law type formulation for 
the medium entropy production in the framework of the 
Onsager-Machlup theory [36, 37], we split the Lagrangian 
C into time reversal symmetric, = C + and anti¬ 
symmetric, = L — C), components and arrive at 




rE = 


1 ■ 
2 . 


- ASm/kB 


(38) 


Using Eqs. (32) and (34), C" in the above equation has 
been replaced by the rate of entropy production ASm- It 
is evident from Eq.(lO) that the average path described 
by 

{x)+ {{x)-vt)/Tr+(j{y) =d, (39) 

(y) + (y) I'Tr -a{x) = 0, (40) 


work done or the energy put into the system by the trap¬ 
ping potential over time [0 : t/] is 

'yVti[x{t)] =—kv f {x — vt)dt. (43) 

Jo 

In order to obtain the distribution of this work, we need 
to perform functional averages of this quantity over all 
possible paths [x{t),y{t)]0 <t<tf] as well as integrals 
over the initial and final points [24, 36, 37] of the path. 
The probability measure of such a path is given by, 

V [x{t), y{t)]0 < t < t/] ~ exp j C{x, y, x, y, t) dt^ . 

(44) 

So, the probability distribution of the work performed by 
the trap, henceforth denoted as Wtr, is 


T’(Wtr) = (^Wtr - J {x- vt) dt 

= ydxoy dyo pixQ.yo) J dxf J dyj 

X / V[x{t)] [ V[y{t)] r[x{t),y{t)] 


X S ^Wtr — fcu J (x — vt) dt 

= ^( dAe*^^'^- P(iAWtr); 

where 

P(iAWtr) = J dxo J dyo p[xo,yo) J dxf J dy/ 

X y 2?[a;(t)] y T»[y(t)] exp J CMix,y,x,y,X)dt^ , 

(46) 


(45) 


corresponds to the minimum value of the integral 
fg-' C{x, y, X, y, t) dt. Substituting (39) and (40) into the 
relation 

i{x), (x), (y), {y),t) = ASm {{x), (x), (y), (y),t) /ks, 

(41) 

we find 

^ [(i) - (( 2 ;) - vt)/Tr - a{y)f 

+ [-{y) + {y)/xr - cF{x)f > 0; (42) 

Eq. (42), ensuring ASm {{■ ■ ■)) > 0, shows the consis¬ 
tency with the second law of thermodynamics. 

5. DISTRIBUTION OF THE WORK 
PERFORMED BY THE MOVING TRAP 

In this section, we derive the distribution of the work 
done by the harmonic trap on the Brownian particle. The 


with 

1 r 2 

CM{x,y,x,y,X) = —t^[x+{x-vt)/Tr + ay] 

+ [ 2 / + y/xr — — 4iAu(a; — ut)/Tr|. (47) 

In order to derive the distribution function at a fi¬ 
nite time, we consider the system to have evolved 
from the initial equilibrium distribution, p(xo,yo) = 

II exp(-^[a;g-kyg]). 

To evaluate the functional integrals in (46), we follow 
the variational approach discussed earlier in section 3. 
We seek the most probable trajectory {x{t), y{t)) which 
corresponds to the largest contribution to the path in¬ 
tegral. The extremization condition, S Cm dt = 0, 
leads to two Euler-Lagrange equations 

X + 2ay — + \It^)x — ^(l — 2iA)/Tr -|- = 0, 

(48) 

y — (cr^ -I- \lTl)y — 2ax + vcrt/Tr = 0, (49) 
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which are supplemented with the boundary conditions, 
x(t = 0) = xq, x{t = tf) = Xf, y{t = 0) = yo and 
y(t = tf) = yf. After some simplifications, the two 
Euler-Lagrange equations can be expressed in terms of 
higher derivatives of a single variable only. These are 

y + 2 (cr^ _ 1 /t^) y + ( 0-2 + Xlrlf y 
— (cr^ + l/r^) vat/Tr + ^avjT^ = 0 , 

(50) 

(cr^ + \ItI) X [y - (cr^ + ^lTl)y + w/Tr] 

+ 2ay — v\/Tr + vt/T"^ , (51) 

where, A = 1 — 2iA. The optimal solutions are 

x{t) =v [bt — Xxr — (A — 2)ct^t(^] 

+ cosh(t/T^) [c 2 cos(tTt) — Cl sin(CTt)] 

+ sinh(t/rr) [c 4 cos(crt) — C 3 sin(crt)] , (52) 


y{t) =aTrv[bt — 2Tr]/b'^ + cosh(t/rr)[ci cos{at) + C 2 sin(CTt)] 
+ sinh(t/Tr)[c 3 cos{at) + C 4 sin(tTt)]. (53) 


The explicit dependence of the constants ci, C 2 , C 3 and 
C 4 on the boundary conditions and other parameters of 
the problem is shown in appendix B. Incorporating the 
quadratic term arising from the fluctuations around the 
most probable path, we write Eq.(46) as. 


P{iXWtT) = j dxo J dyo p{xo,yo) j dxf 


X / dyj exp[-J CM{S:,y,x,y, X) dt ] . 


(54) 


Substituting the solutions (52) and (53) into (54) and 
subsequently performing first the integral over time and 
then over final and initial positions (some of the interme¬ 
diate calculations are shown in appendix B), we have 

P{iXWtr) = exp [- (VEoA^ -f ilTA)] , (55) 

where 

+ cos{(Ttf) \Tr — — 2crr2 sin(crt/)|; (56) 

W = ^ 2Db^ [-2.Tr + 6cr^r)^ -f a’^b’^tjvr 

+ 2 tf(l — cr^T^)] + cos{atf) [ 2 Tr — 

-|-2crr^ sin(crt/) [cr^r^ — 3] (57) 


Substituting (55) into Eq.(45) and doing the final Gaus¬ 
sian integral over A, we find the final form of the work 
distribution as 


1/2 


exp 


iPWtr - tT)2 


4Wo 


• (58) 



FIG. 2: Plot of the symmetry function with Wtr for the 
parameter values D = 0.01, /3 = 100, Tr = 0.5, k = 2, 
and V = 2, aX the final time tf = d. The dashed-dot line 
corresponds to the situation where the symmetry 
function follows the fluctuation theorem. The deviation 
from the fluctuation relation increases for larger values 

of a. 


The known result of the work distribution function for 
a uniformly dragged colloidal particle [18, 19, 21-23] can 
be retrieved by substituting ct = 0 in our general result. 
The work distribution function, in this case, obeys the 
fluctuation theorem. The symmetry function defined as 


Sw 


PjWtr 
P(-Wtr)J ’ 


therefore, has a value Sw 


PWtr. 


In the presence of the non conserving force (cr ^ 0), the 
symmetry function is 


Sw = pWtrWIWo. (59) 


Clearly it does not satisfy the conventional fluctuation 
theorem since to hold this relation, it is required that 
Sw = fdWtr- The present work distribution function sat¬ 
isfies the fluctuation theorem if Wq = W. Given the 
following expressions for Wo and W 


V 


w ^ 


(60) 


2Db^ + 6cr^T(? -I- a'^b'^t'fTr -f 2t/(l - , 

(61) 


in the steady state limit tf/xr >> 1 , such a condition is 
satisfied if 

a'^Tr [6t^ — 2 cr^r(^ + t'jb'^ — dtfTrb] = 0. (62) 

tr = 0 trivially satisfies this equation. The other nonzero 


P(VEtr) 


47rWo 


















values of a for which Eq. (62) is satisfied are 


2tfTr -t}+ ^(- 4(^)2 + 4 ^ + 1 ) 

(63) 

1 

2tfTr -t)+ t 2 + ,/(- 4(^)2 + 4 ^ + 1 ) 

(64) 

In the steady state limit {tfjTr » 1), these roots be¬ 
come complex and hence not physically acceptable. To 
see the deviation from the fluctuation theorem, we have 
plotted the symmetry function as a function of Wtr for 
different values of a (see figure (2)). The deviation be¬ 
comes more prominent as one increases the value of a. 


1 

CTl = -- 

IfTr 

= -cr2; 

1 

= - - 

IfTr 

= - 0 - 4 . 


6. SUMMARY AND PERSPECTIVES 

In this article, we have studied the two dimensional 
motion of a Brownian particle subject to a time depen¬ 
dent harmonic force and a nonconservative rotational 
torque. Starting with the overdamped Langevin equa¬ 
tions, we have obtained a path integral description for 
the transition probability of the Brownian particle. The 
transition probability gives the probability of finding the 
Brownian particle at a final position at time t / given that 
it was at a specific position at an earlier time to- Us¬ 
ing a variational method, we evaluate the path integral 
to obtain the explicit form of the transition probability. 
This, upon integration over the initial position coordi¬ 
nates gives the joint probability distribution. 

In order to obtain the definitions of the work, the dis¬ 
sipated heat and the entropy production of the medium, 
we consider the forward and time-reversed trajectories of 
the particle. Here, the time reversed trajectory is char¬ 
acterized by the reversal of the velocity variables and 
the flipping of the end points with the same underlying 
dynamics. After identifying the medium entropy pro¬ 
duction and the dissipated heat, we obtain the energy 
conservation relation and, subsequently, define the two 
works; one done by the moving trap and the other by the 
nonconservative force. Further, we show that the aver¬ 
age path corresponding to the minimum of the general¬ 
ized Lagrangian leads to the positive entropy production 
rate. 

The expression of the thermodynamic work further al¬ 
lows us to derive its distribution function using a similar 
path integral approach adopted for the transition proba¬ 
bility. We expressed the distribution function in a func¬ 
tional integral form which was approximated by the max¬ 
imum contribution coming from the optimal path and the 
fluctuations upto second order about the optimal path. 
The distribution is Gaussian in nature. The explicit form 
of the distribution correctly produces the known result of 
the work distribution for a dragged Brownian particle in 
the absence of any other force (cr = 0). However, in the 


presence of the nonconservative force, its mean and vari¬ 
ance have a non-trivial oscillatory dependence with a as 
the angular frequency. As a consequence, the work distri¬ 
bution fails to satisfy the conventional work fluctuation 
relation. 

In this context, we discuss another version of the model 
where the thermodynamic work naturally appears as a 
ratio of the forward and reverse process and satisfies fluc¬ 
tuation relation. Here, the nonconservative force is ap¬ 
plied at the position of the minimum of the moving har¬ 
monic trap. In the reference frame of the trap, the parti¬ 
cle moving under the potential cj}{X,y) = k{X^ + y^)j2, 
satisfies the overdamped Langevin equations given by, 

X = -Xj Tr-v - ay + ix/oi, (65) 

y = -yj Tr^aX + ^y/a, (66) 


with, X = X — vt. 

To obtain the thermodynamic work and its fluctuation 
relation, we consider the time reversed process with the 
reversal of sign in appropriate terms of the above equa¬ 
tions [41]. The corresponding Lagrangian, denoted by 
L^, has the form 




1 

Id 


{-X + X/Tr -V- ayf + {-y + yj 


Tr 


+ aXf 
(67) 


It is straightforward to obtain the ratio of the forward 
and reverse transition probabilities as 


El 

Pr 


= exp 


P 


(^—X(j) — kv 



( 68 ) 


The first term in the argument of the exponential stands 
for the change in the potential energy and the second 
term can be recognized as the thermal work done [8]. If 
we consider the initial distribution to be a canonical one, 
i.e. po oc then the ratio of forward and re¬ 

verse probability densities yields the fluctuation relation 
for the work >Vtr = / —kvX dt, 

^ = exp (/JVVtr) , (69) 

Pr 

where the probability measure of the path starting at a 
given initial point is obtained from the transition prob¬ 
ability as 'Pf,r = Po Pf,r- In this context, it will be of 
interest to verify whether the distribution of this work 
satisfies various aspects of fluctuation relations and will 
be discussed elsewhere. 

Finally, the model studied in this paper and its exten¬ 
sion seem to be promising for further investigations and 
the results obtained would also be amenable to experi¬ 
mental verifications. 
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Appendix A: Path integral solution of the Fokker-Planck equation 


The joint probability density function p[x, y, < + e) at time t + e is related to the transition probability 
P{x, y,t + e \ x', y', t) as 


p(x, y,t + e) = J P{x, y,t + e\x', y', t) p{x', y', t) dx' dy'. 


(Al) 


For infinitesimal time difference e, the transition probability can be expressed as, 


P{x,y,t + e\x',y',t)^ 1 + Lfp • e + O(e^) 6{x - x') S{y - y') 


(A2) 


Using the Fourier representation of the S function and substituting App from Eq.( 6 ), we have, 


^9 ^52 


P{x,y,t + e\x',y',t) = J dqi J dq 2^^1 + e{{x'- vt)/Tr + ay') 

d 52 'I 

-£{-y'lTr + j exp[i(gi(a: - a:') +g 2 (y - y'))] 


-AD 


(x — x') 


+ {x' — vt) /r-T + ay' 


e + 


iy-y') , 


+ y'/Tr - ax' 


(A3) 


In order to obtain the transition probability for a finite time difference (t/ — to), we divide the time difference into 
N small segments of duration e = {tj — tQ)/N . Then, in the limit A —>■ 00 and e ^ 0, the repetitive application of 
the Chapman-Kolmogorov equation [44] yields, 


f ^ dx- f ^ 

p(x,,„,o I = J n J n ,4,^,,,./. 


dy. 


f ^ 1 1 

X expi - X! ^ “ Xi-i)/e + (xj - vti)fTr + cry*]^ + [(y* - yi-i)/e + yt/xr - axif^ e I 

J '^[xit)] J P[y{t)] expj- ^ ^ (^[x + (x - vt)/Tr + ayf + [y + y/r^ - axf^ dtj. 


(A4) 


Appendix B: Some intermediate calculations of section 5 


After substituting the optimal solutions (52) and (53) into Eq.(47), the expression of the modified Lagrangian reads 

-|^ 2^2 + X)vTr - e‘/'^’’ 5 [(c 2 + C 4 ) cos(tTt) - (ci + C 3 ) Sm((Tt)] 

r{(-l + A)tTr2z; + e ‘/'^’’6 [(ci + C 3 ) cos{at) + (c 2 + C 4 ) sin(CTt)]}^ 

+ " „ ^ i —Trv[Xb + a^txrb — 2a‘^T^]l\?' + cosh(t/rj.)[c 2 cos{at) — a sin(crt)] 


£M(i,y,i,y, A) = 
+ 


dDr^b^ 

1 

4iAr262 

(A - 1 )^; 

2DTr 


+ sinh(t/T4.)[c4 cos((Tt) — C3 sin(crt)] >. 


(Bl) 
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where, Ci, C2, C3 and C4 have the form 

Cl =yo + ^aT^vjb'^, 

C2 = Xo + VTr{Xb — 2(P‘Tl)lh^, 

C 3 = — coth(t//rr)[yQ 6 ^ + 2aT‘^v] + cschitf/Tr) '[—aTrvitfb — 2Tr) cos{atf) + yfb^ cos(crt/) 

+ V sm{atf){tfb — Xrrb + 2a^T^) — Xfb^ sin(crt/)] 

C4 = ^ 1 — coth(<//rr)[t'Tr(A6 — 2 a^T^) + x^b^] + csch(t//Tr) [vTr{Xb — 2 a^T^) cos{atf) + {xfb^ — vtfb) cos{atf) 
+ yfb^ sin(CTt/) — aTrv{tfb — 2tj.) sin(CTt/)] 


(B2) 

(B3) 


(B4) 


(B5) 


We use the above expressions in Eq.(54) to evaluate various integrations. After performing the time integration 
and integrations over the final positions, we have 

P(iAWtr) = J dxo J dyo p{xo,yo) - X)[hiXo + h2yo + hsl/idDb^)} 


exp(^—ve — A)/i3/(4Z)6^)^ J da;oexp|— pkxl/2 + ve — X)hiXQ/{ADb^) 


27r 

X / dyo exp<j - 


Pkyl/2 + ue-‘^/"’'(l - X)h2yo/idDb^) 

where /ii, /12 and ho are 

hi = — 26 ^ [cosh{tf/Tr) — cos{atf) + ar^ sin(crt/) + sinh(tj/rj.)], 


h2 = -2b^ 


(jTr cos{atf) — + sin{atf) 


and 


(B6) 

(B7) 

(B8) 


ho = 2Tf.ucos((Tt/)[A6 — 4 (t^t^] + v cosh(t//rr) [Scr^r^ — 2XTrb + a^tj^Trb^ + btf{l + bX — Scr^r^)] 

— 2(Tr^usin((Tt/)[2 + bX — 2a^T^] + usinh(t//rr)[ 6 cr^r^ ~ (1 + X)Tr + (A — 

+ + t/5(l + A 6 — 3(t^t^)] . (B9) 

After doing the Gaussian integrals in Eq.(B 6 ) and using the relation l/{D‘^j3k) = Tr/H, we arrive at Eq.(55) in the 
main text. ,_ 
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